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Abstract 

In this paper, by using bifurcation method, we successfully find the Fornberg- 
Whitham equation 

V>t U xx i -\- U x — WUxxx ~\~ 3u x U XX i 

has a type of traveling wave solutions called kink-like wave solutions and antikink- 
like wave solutions. They are defined on some semifinal bounded domains and pos- 
sess properties of kink waves and anti-kink waves. Their implicit expressions are 
obtained. For some concrete data, the graphs of the implicit functions are displayed, 
and the numerical simulation is made. The results show that our theoretical analysis 
agrees with the numerical simulation. 
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1 Introduction 



It is well known that the exact solutions for the nonlinear partial differential 
equations can help people know deeply the described process. So an important 
issue of the nonlinear partial differential equations is to find their new exact 
solutions. Traveling wave solution is an important type of solution for the 
nonlinear partial differential equations and many nonlinear partial differential 
equations have been found to have a variety of traveling wave solutions. For 
instances, the well-known Korteweg-de Vries equation 

u t — Quu x + uu xxx = (1-1) 

has solitary wave solutions and its solitary waves are solitons [1]. A KdV-like 
equation 

u t + a(l + bu n )u n u x + 5u xxx = (1.2) 

has some kink wave solutions [2]. Its kink wave solution = x — ct) was 

difined on (— oo, +oo), and lim u(£) = A, lim = B, where A and B are 

oo £— >oo 

two constants and A ^ B. The Camassa-Holm equation 

Uf — ^xxt ^>uu x — *2u x u xx ~\- uu xxx (1.3) 

has peakons, cuspons, stumpons, composite wave solutions [3], it also has 
compactons [4]. The Degasperis-Procesi equation 

u t - u xxt + Auu x = 3u x u xx + uu xxx (1.4) 
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has a multitude of peculiar wave solutions: peakons, cuspons, composite waves, 
and stumpons [5]. The Kuramoto-Sivashinsky equation 

u t + uu x + u xx + u xxxx = (1.5) 

has periodic and solitary solutions [6]. In [7], Liu, Li and Lin found a new type 
of traveling wave solutions for the Camassa-Holm equation, which are defined 
on some semifinal bounded domains and possess properties of kink waves or 
anti-kink waves. They called them kink-like waves and antikink-like waves. 
Later, Guo and Liu [8] found the CH— 7 equation 

Ut + c u x + 3uu x - a 2 (u xxt + uu 

xxx 3u x U xx ) -\- r )U xxx 0, (1.6) 

posses kink-like waves when a 2 > 0. Tang and Zhang [9] showed the CH— 7 
equation (1.6) also has kink-like wave solutions even when a 2 < 0. Chen and 
Tang [10] showed that the Degasperis-Procesi equation (1.4) has such type 
of traveling wave solutions. Recently, Liu and Yao [11] found the following 
generalized Camassa-Holm equation 

Ut + 2ku x — u xxt + auu x = 2u x u xx + uu xxx (1-7) 

also posses kink-like wave solutions. 

We are motivated to seek kink-like wave and antikink-like wave solutions 
for the Fornberg-Whitham equation 

Ut — u xxt + u x = uu xxx — uu x + 3u x u xx . (1-8) 

To our knowledge, such type of traveling wave solution has never been found 
for the Fornberg-Whitham equation. Eq.(1.8) was used to study the qualita- 
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tive behaviours of wave-breaking [12]. It admits a wave of greatest height, as a 
peaked limiting form of the traveling wave solution [13], u(x, t) = A exp(- 
where A is an arbitrary constant. 



x 3 t 



The remainder of the paper is organized as follows. In Section 2, we state 
the main results which are implicit expressions of the kink-like wave and the 
antikink-like wave solutions. In Section 3, we give the proof of the main re- 
sults. In Section 4, we make the numerical simulation of the kink-like and the 
antikink-like waves. A short conclusion is given in Section 5. 

2 Main results 

We state our main result as follows. 
Theorem 1. For given constant c, let 

i = x-ct (2.1) 



^ = c-l±^(c-iy-2g (2.2) 
9i(c) = (2.3) 
*<c) = (2-4) 



(1) If g < (72(c), then Eq.(1.8) has two kink-like wave solutions u = y?i(£) and 
u = ^3(0 an d two antikink-like wave solutions u = <^ 2 (0 an d u = ^(O- 



{2^ipl + l 1 ip l + l 2 + 2i Pl + l l )( Vl -ip^ _ L 

= p ie »«, £e (-00,$) (2.5) 



{2^/ElJipl + l 1(pi + l 2 + b m + Z 3 ) Ql 



— \ =/3ie» 6 , £ £(-£„, 00) (2.6) 

(2y/aiJ<pl + h<P2 + k + bnp 2 + / 3 ) Q1 
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— \ = = /3 2 e » e ,£e(-$,oo) 

(2 v / a^^/( / 93 + miip 3 + m 2 + 62^2 + m 3 ) a2 

(2.7) 



and 



(2 Jv9| + miV9 4 + m 2 + 2ip 4 + mi)(^ 4 - V^o)" 2 „ i f _ - 

= /? 2 e^,e e (-00, e 



where 



/ 2 = - [1 - 4c + 3c 2 - 3<? + (3c + 1) ^(c - l) 2 - 2g] 



h = 3 [2 - 5c + 3c 2 - 6,7 + (3c + 2) ^(c - l) 2 - 2<j] 



m 1 = -(l-3c + 3 v /(c-l) 2 -2( ? ) 

m 2 = -[1 - 4c + 3c 2 - 30 - (3c + 1) yf(c - l) 2 - 2g] 



m 3 = - [2 - 5c + 3c 2 - 6g - (3c + 2) yj(c - l) 2 - 2</] 



at = 4(1 - 2c + c 2 - 2# + J(c - l) 2 - 2g) 



a 2 = 4(1 - 2c + c 2 - 2g - J(c - l) 2 - 2^) 
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6 1 = ---4 v /(c-l) 2 -2 5 



i 



6 2 = -- + 4 v /(c-l) 2 -2( ? 



1 + ^(c-l) 2 -2^ 

ai = — 



2 J( c -l)*-2g + J{c-l)*-2g 



(2.8) 



/ 1 = -(l-3c-3 v /(c-l) 2 -2( 7 ) (2.9) 



2.10) 
2.11) 
2.12) 
2.13) 
2.14) 
2.15) 
2.16) 
2.17) 
2.18) 

2.19) 
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:-l) 2 -2o 

a 2 = vv — (2.20) 



2^{c-lf-2g-^c-lf-2g 
go = (2V^ + hc + l 2 + 2c + h)(c-^ 

Pl (2^Vc 2 + hc + l2 + b lC + l 3 )^ { ' ' 

o = (2 Vc' + miC + W 2 + 2c + m 1 )(c- <p^T (2 22) 

2 (2^Vc 2 + rrhc + m 2 + b 2 c + m 3 ) Q 2 

_ ^ (2y^ + 7^ + 7; + 2a + ji) (a - ^ 
2y/aiy/a 2 + ha + l 2 + b x a + Z 3 )<w 

= (2y/fe 2 + m x & + ^~ 2 + 2b + mi) (6 - 24 
2 {2y/aiy/V l + m x b + m 2 + b 2 b + m 3 )«2 

e 1 = 21n(/5 1 //5 1 °) (2.25) 
£ 3 = 21n(/3 2 °//3 2 ) (2.26) 



a and b are two constants satisfying <fii(0) = y? 2 (0) = a, ^(O) = 924(0) = b, 
and (/9g<a<c<6< • 

(2) If g 2 (c) < g < g\(c), then Eq.(1.8) has a kink-like wave solution u = <^i(0 
of implicit form (2.5) and a antikink-like wave solution u = ip 2 (Q °f implicit 
form (2.6). 

We will give the proof of this theorem in Section 3. Now we take a set of 
data and employ Maple to display the graphs of u — <fi(£)(i = 1, 2, 3, 4). 

Example 1. Taking c = 1 and g = —4 < g 2 (c) (corresponding to (1) of 
Theorem 1), it follows that </?o = -2.82843, tp$ = 2.82843, h = -6.99019 and 
l 2 = 15.5425, h = 50.8562, a x = 43.3137, h = -12.647 , a x = -0.581712. 
Further, choosing a = —1 G (<Po,c), we obtain £q = 0.387475. We present 
the graphs of the solutions <pi(£) and <f 2 (C) hi Fig-1 ( a ) an d (b) respectively. 
Meanwhile, we get mi = 4.32352 and m 2 = 0.457528, m 3 = 13.1438, a 2 = 
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1 2 3 4 5 6 -6 -5 -4 -3 -2 -1 



6 5 
(c) (d) 
Fig. 1. The graphs of <pi{£)(i = 1, 2, 3, 4) when c = 1, g = -4, a = -1, b = 2. 
20.6863, b 2 = 9.98038, a 2 = 0.40201. Further, choosing b = 2 G (c,^), we 
get £q = 0.274787. The graphs of the solutions ^(O and ^(O are presented 
in Fig. 1(c) and (d) respentively. The graphs in Fig.l show that <fi(£) and 
</53(0 are kink-like waves and ^(O and 934 (£) are antikink-like waves. 

3 Proof of main results 

Let u = with £ = x — ct be the solution for Eq.(1.8), then it follows 
that 

- ap' + op'" + tff = pp'" - pp' + 3p'p" (3.1) 

Integrating (3.1) once we have 
p"(p - c )=g-cp + p + ^p 2 - {p'f (3.2) 
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where g is the integral constant. 



Let y = ip', then we get the following planar system 



^ = v 



d£ ip—c 



(3.3) 



with a first integral 

Hfay) = (<p- c)V - - \{<p - c) - g - c+ ^] = h (3.4) 



where h is a constant. 

Note that (3.3) has a singular line ip = c. To avoid the line temporarily 
we make transformation d£ = (ip — c)d(. Under this transformation, Eq.(3.3) 
becomes 

% = (<P-C)V 

(3.5) 

f^g-ap + V + ^-V 2 

The Eq.(3.3) and Eq.(3.5) have the same first integral as (3.4). Consequently, 
system (3.5) has the same topological phase portraits as system (3.3) except 
for the straight line tp — c. Obviously, <p = c is an invariant straight-line 
solution for system (3.5). 

Now we consider the singular points of system (3.5) and their properties . 
Note that for a fixed h, (3.4) determines a set of invariant curves of (3.5). As 
h is varied (3.4) determines different families of orbits of (3.5) having different 
dynamical behaviors. Let M(tp e ,y e ) be the coefficient matrix of the linearized 
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system of (3.5) at the equilibrium point (<p e ,y e ), then 
/ \ 

He <£e-C 



M(p e ,y e 



ip e -(c- 1) -2y e 



and at this equilibrium point, we have 

J(ip e ,y e ) = det M(p e ,y e ) = -2y 2 - (p e - c)[ip e - (c - 1)], 
p(p e ,ye) = trace(M(<p e ,y e )) = -y e . 

By the theory of planar dynamical system (see [14]), for an equilibrium point 
of a planar dynamical system, if J < 0, then this equilibrium point is a saddle 
point; it is a center point if J > and p = 0; if J = and the Poincare index 
of the equilibrium point is 0, then it is a cusp. 



Since system (3.3) has the same topological phase portraits as system (3.5) 
except for the straight line p = c. By investigating the topological dynamics 
of system (3.5), we can obtain the following properties for system (3.3). 

(1) If g < (72(c), then system (3.4) has two equilibrium points (y?o,0) and 
(<Pq , 0). They are saddle points and there is inequality </?q<c — 1 < c < <Pq. 
In this case, there are four orbits connecting with ((^g,0), we use and 
I 2 - to denote the two orbits lying on the right side of ((pn,0) (see Fig.2(a)). 
Meanwhile, there are four orbits connecting with (y9g",0), we employ /*+ and 
l 2 + to denote the two orbits lying on the left side of {(fo, 0) (see Fig. 2 (a)). 

(2) If (72(c) < g < <7i(c), then system (3.4) has two equilibrium points (ip^ , 0), 
(<Pq, 0). (<Pq , 0) is a saddle point and {ip^ , 0) is a center point or a degenerate 
center point. and y?d satisfy that (p^ < c — 1 < ip^ < c. and l 2 - are 
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used to denote the two orbits lying on the right side of (o9 , 0) (see Fig.2(b)) 
(3) If gi(c) < g, then system (3.4) has no equilibrium point. 



1 


<p = C 




l\ 














/ 2X 


/ 1' 



<p 





1 <p = 
















» 


/ 2X 





<Po <Po 

(a) g< 52(c) (b) 52(c) < g < 51(c) 

Fig. 2. The sketches of orbits connecting with saddle points. 



On the 09 — y plane, the orbits I _, I _, I + and I + have the following 
expressions respectively. 



,1_ = 1 (V?-<A))vV 2 + ^ + / 2 , 36 ^ 

' y 2 c - 09 1 ' ' 



/2 _ = 1 (v?o -y)vV 2 + /iy? + / 2 

^0 ' y 2 c - 09 1 ' j 

,1 1 (v?o ~ V?)vV + mnp + m 2 

' y 2 09 - c 1 ' ' 

/2 1 {<P ~ <foW^_ + mnp + m 2 

' V 2 09 - c 1 ' J 



where o9 and o9q~ are in (2.2), l\ and / 2 are in (2.9) and (2.10), m\ and m 2 are 
in (2.12) and (2.13). 

Assume that tp = (p = tpn(£), ip = tp 3 (£) and 09 = (p 4 (£) on I 1 I 2 

/* + and I 2 + respectively and o9i(0) = y9 2 (0) = a, oj 3 (0) = o9 4 (0) = b, where a 
and 6 are two constants, and a £ (o9q ,c), b £ (c, o9g). Substituting (3.6)-(3.9) 
into the first equation of (3.3) and integrating along the corresponding orbits 
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respectively, we have 

/- N C ~ S =ds = - [ ds (along I 1 -) (3.10) 









/- N C y * =t/s = -[ ds (along Z 2 _) (3.11' 

b - 1 ° 

fr+ — rkr r ds = -Jds (along i\a (3.12) 

^=rfg = - / ds (along Z 2 + ) (3.13) 

(s — LPn )\/s z + mi s + mo 2 J v o 







Computing the above four integrals we obtain the implicit expressions of </?i(0 
as (2.5)-(2.8). 

Meanwhile, suppose that y?i(£) — > c as £ — > £q, ^(O — > c as £ — > — £ 2 , 
^3(0 — > c as £ — > — £q, (^4(0 ~^ c as £ ~^ £0' then it follow from 
(3.10) -(3.13) that 

c 

£0 = £ 2 = / 7 — M , ds ( alon s Co) ( 3 - 14 ) 

J (S - Lf )y/S 2 + IjS + Z 2 ^° 



^ = ^0 = / 7 + , T = ds ( al0n S ( 3 - 15 ) 



b 



Computing the above two integrals, we get the expressions of £q and £q as in 
(2.25) and (2.26). The proof is finished. 



4 Numerical simulations 



In this section, we will simulate the planar graphs of the kink-like and the 
antikink-like waves. 
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From Section 3, we see that in the parameter expressions (f = and 
V = 2/(0 °f the orbits of System (3.3), the graph of and the integral 
curve of Eq.(3.2) are the same. In other words, the integral curves of Eq.(3.2) 
are the planar graphs of the traveling waves of Eq.(1.8). Therefore, we can 
see the planar graphs of the kink-like and the antikink-like waves through the 
simulation of the integral curves of Eq.(3.2). 

Example 2. Take the same data as Example 1, that is c = 1, g = —4, 
a = — 1, b = 2. Let (p = a = — 1 in (3.6) and (3.7), then we can get y ~ 
2.217532085 or y « -2.217532085. And let <p = b = 2 in (3.8) and (3.9), then 
we obtain y « 1.499467912 or y « -1.499467912. Thus we take the initial 
conditions of Eq.(3.2) as follows: (i) Corresponding to I 1 - we take (p(0) = 
-1 and <ff(0) = 2.217532085. (ii) Corresponding to Z 2 _ we take <p(0) = -1 
and ip'(0) = —2.217532085. (iii) Corresponding to /*+ we take y?(0) = 2 and 
ip'(0) = 1.499467912. (iv) Corresponding to l 2 + we take ip(0) = 2 and ip'(0) = 
-1.499467912. 

Under each set of initial conditions we use Maple to simulate the integrals 
curve of Eq.(3.2) as Fig. 3. Comparing Fig.l with Fig. 3, we can see that the 
graphs oitpi(£)(i = 1, 2, 3, 4) are the same as the simulation of integrals curve 
of Eq.(3.2). This implies that our theoretic results agree with the numerical 
simulations. 

5 Conclusion 

In this paper, we find new bounded waves for the Fornberg-Whitham equa- 
tion (1.8). Their implicit expressions are obtained in (2.5)-(2.8). From the 
graphs (see Fig.l) of the implicit functions and the numerical simulations (see 
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(c) ip(0) = 2, p'(0) = 1.499467912. 



-3 
xi 




(d) ^(0) = 2, p'(0) 



-1.499467912. 



Fig. 3. The numerical simulations of integral curves of Eq.(3.2). 
Fig. 3) we see that these new bounded solutions are defined on some semifinal 
bounded domains and possess properties of kink waves and anti-kink waves. 
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